Exact solutions in Einstein cosmology with a scalar field 
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We present a simple way to obtain exact solutions of Einstein-scalar field equations on spatially 
flat Friedmann-Robertson- Walker space-times. The scalar equation turns out to be integrable if 
the Hubble parameter is written as an appropriate function of the scalar field and its velocity. 
Eventually, the field equations are reduced to find 'generating functions' for a given scalar potential. 
Once a generating function is found as a function of the scalar field, the evolution of the field and 
the Universe can be easily obtained with a simple integration. As examples, we obtain the solution 
spectra in the cases of the constant and the exponential potentials, and find exact solutions for 
various scalar potentials such as the Xcj) 4 , the power law, and the double-well hyperbolic functions. 
We additionally analyze the stability of the generating equation. We show that the existence of 
a fixed point of the equation of motion affect on the evolution so that the Universe experiences a 
long inflation. We additionally show that small change of the scalar potential cannot get rid of the 
appearance of the long inflation. 
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INTRODUCTION 



O h Historically, the scalar field is one of a main elements in cosmology used as the inflaton [l| , as the dark matter 0] , 
and as the dark energy candidate J3( • The scalar field coupled to gravity have also appeared in connection to the so 
bJ[) called string theory landscape [3, [5j], where the scalar potential V(4>) is usually thought of as having many valleys, 
which represent the different vacua solutions. The dimensional reduction of higher dimensional fundamental theory to 
^ four dimensions typically gives rise to scalar fields with exponential potentials coupled to four-dimensional gravity . 
Recently ATLAS and CMS at the Large Hardron Collider reported a discovery of Higgs-like particle p}, which 
motivates researches on the scalar field. If it is a Higgs in the standard model then the mass and quartic coupling 
■ of the Higgs potential are now determined with all the interactions of Higgs to the other standard model particles. 
The Higgs boson was discussed as a seed of inflation Q , in which the scalar potential has a broken symmetry with 
</> 4 interactions. Even though it was argued that there are unitarity problem and the instability of the potential up 
to the near Planck scale @, the Higgs may have played a major role to seed the formation of the structures in the 
present unverse [lOj . 

In this work, we are mainly interested in finding exact cosmological solutions of Einstein equation coupled to a 
scalar field <fi with action in standard form, 



(1) 



where we set Mpi = 1, h = 1 = c. The universe is spatially flat, homogeneous, and isotropic, with metric: 

ds 2 = -dt 2 + a 2 (t)(dx 2 + dy 2 + dz 2 ), (2) 
where a(t) is the scale factor. The dynamics of the scalar field and gravity can be dealt with a pair of equations 



8ttG 



+ 3fty + V"($ =0, (4) 

where dot and prime denote derivatives with respect to time and the scalar field, respectively, and H — a/a is the 
Hubble parameter. The time derivative of Eq. ([3]), by using Eq. leads to a kind of Riccati equation for the Hubble 
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parameter, 3H 2 + H = V((f>). Methods relating the one dimensional time-independent Schrddingcr equation to the 
Riccati equation were developed [Til ]. In the case of an exponential potentials, the scalar cosmology in four dimensions 
were investigated [l2| and general exact solutions were found [l3[ . Interesting properties of the cosmological solutions 
with exponential solutions were also discussed fl4l - [l6j . In a (phantom) scalar-tensor theory, the late-time cosmology 
was studied with an exponential potential and by using reconstruction technique [l6j . For the case of tachyon scalar 
scalar field, Padmanabhan (l7j have shown that one can reconstruct a corresponding potential once a time-dependent 
scale factor is given, which result could be extended to general cases. 

In this work, we present a new simple way to find exact solutions of the differential Eqs. ([3]) and (HJ). In Sec. II, 
the equation of motion is reduced to a 'generating equation': a problem finding the generating function through a 
non-linear first order differential equation. The evolution of the scalar field and the scale factor are related with the 
generating function in a simple manner. In Sec. Ill, we obtain the generating function for the cases of the constant 
and the exponential potentials by attacking the generating equation directly and find the evolutions of the universe. 
However, in general, the generating equation is hard to solve and we detour the difficulty in Sec. IV by choosing the 
generating function first and obtain the scalar potential algebraically later. As examples we present exact solutions 
for several different potentials. We additionally analyzed the properties of the fixed point of the generating equation 
and show that its property gives a long inflating period on the Universe. In Sec. V, we show that the long inflation 
is a frequent occurrence in the evolution of the Universe by showing that it happens even the initial conditions or the 
potential are perturbed. In Sec. VI, we summarize and discuss the results in relation to the inflation. 



II. THE GENERATING FUNCTION FOR THE SCALAR COSMOLOGY 



The system of the coupled equations ([3]) and ((4]) has two unknowns, cf>(t) and H(t) for a given potential V(4>). It 
is easy to see that the equation Q is integrable if we set, 

where G{<f>) is an arbitrary function of the field, which we call 'generating function'. Now, the scalar equation of 
motion (|4|) is integrated to give T = \(\? — G 2 {<f)) ~ V((f>), where an integration constant is absorbed into the 
definition of G(4>). Using this result, the Einstein equation ([3]) becomes the 'generating equation': 

U(0) = G 2 (0)-^[GW, (6) 

where we have removed the trivial solution G = leading the fiat space-time. Using ©, the scalar field evolution 
equation and the Hubble parameter are given by 

*-£-><*>. m 

To have an expanding Universe, we should choose the sign of G{4>) to be non-negative 1 . 

As a result, the two coupled differential equations (PJ) and (01 with respect to time is reduced to one non-linear first 
order differential equation ([5]) with respect to the scalar field supplemented by the equation giving the dynamics (J7])- 
If we solve Eq. ^ for a given potential V(4>) and obtain the 'generating function' G(0), the whole solution spectra 
can be found. We succeed in solving Eq. (|6|) for the cases of the constant and the exponential potentials. For most 
cases other than the two, Eq. © is too hard to attack directly. Therefore, we detour the difficulty by specifying the 
generating function first and determine the potential algebraically from Eq. The time evolutions of the scalar field 
and the Hubble parameter are simply given from Eq. (J7]). Since G((j>) and G'(cj)) are directly related to the Hubble 
parameter and the scalar field evolutions, we may control the characteristic behaviors of the solutions easily. As a 
return, we cannot control the form of the scalar potential completely 

Before end this section, we display the acceleration of the scale factor and the equation of state parameter of the 
scalar field during the evolution in terms of the generating function. The acceleration is given by 

-=ij + i? 2 = ^[G(0) 2 -2G'(0) 2 ]. 



1 In Ref. Ilgjl , Reyes found an algebraic way to find the solutions of the scalar cosmology which is close to us. In their work, he found 
several solutions corresponding to power, hyperbolic, and Morse type potentials. From private communications, we notice that similar 
ways to Eq. (£6j) is called as a fake supergravity method [l9l42ll . 
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Therefore, for the region where G(4>) 2 > 2G'(4>) 2 the universe will be expand with accelerating rate. The equation of 
state parameter of the scalar field becomes 

p 4GW 

At the fixed point of Eq. © satisfying G(<fi) 2 = V(4>), the equation of state becomes w = — 1 and the scalar field will 
behaves as if it is a cosmological constant. 

III. EXACT SOLUTIONS FOR CONSTANT AND EXPONENTIAL POTENTIALS 

In this section, we directly solve Eq. ([6]) to obtain the generating function in the cases of the constant and exponential 
potentials. We additionally present the exact solution spectra for the two potentials. 

A. Constant potential 

We first consider the simplest constant potential, 

V(4>) = A > 0. 

For the constant potential, Eq. ((6j) has two different generating functions as solutions. The first solution is G{4>) = VA 
which is the fixed point solution and makes the scalar field does the same role as the cosmological constant. Then, 
we have </> = and H = Hj = y/A/3. This gives nothing but the de Sitter space-time with its scale factor expanding 
exponentially. 

The second generating function is given by 

GW = ~ • 



The scalar field and the scale factor behave as 

RTTr . \ 

(8) 



y|log ^V3iJ/tanh(^-t)^ , a(t) = a sinh 1/3 (3 J ff/t). 



The domain of time is (0,oo). The scalar field monotonically decreases from infinity at t = to a constant value, 
g log A. The universe approaches the de Sitter state of the fixed point. 

B. Exponential potential 

We next consider the exponential potential, 

V{4>) = hs^M. (9) 

Without loss of generality we set (3 > 0. Since the derivative of the exponential is nothing but an exponential, one 
may easily guess a generating function of the form: 



^) = y r ^ev / f^, (10) 

where a real generating function of this form exists only when \/3\ < 1. The scalar field and scale factor corresponding 
to this is given by 



a(t) = a Q (l + p*J-J±^t 



(ID 
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The whole cosmological solutions with the potential ((9j) were studied in Ref. [22j by changing the equations into two 
Riccati equations using a couple of coordinates transformation. The model was also studied in terms of Noether charge 
method in Ref. [23| and Hamilton- Jacobi method ■ The model is extended to include a perfect fluid numerically in 
Ref. 0J25I]. The solution pT|) corresponds to a solution approaching to the fixed point (x, y) = (X/y/6, (1 — A 2 /6) 1 / 2 ) 
in Ref. [1J]. 2 Another two fixed points (±1,0) in the reference corresponds to V = 0, which is not relevant at the 
present situation. 

In this work, we present the whole generating functions for the exponential potential by solving Eq. (|6|) directly. 
The general solution of Eq. ([5]) representing expanding universe is 



V-W = GW = ^ Y -± m eVi^ (12) 

where <p c is a parameter characterizing the maximum possible value of the scalar field during evolution, which will be 
shown below, and the function L(<p) is implicitly defined by the relation 

ll-il^ll + il^- 1 ' { ' 

For j3 = 1, this equation is ill-defined. But, by using (3 — > 1 limit after taking the logarithm at both sides, we get 

l+L ' y ' 

Therefore, we do not deal this case separately. The solution pT|) corresponds to the limit <f> c — > 00 with L((f>) = /3. 
The evolution of the scalar field is given by 

where we have used dL J£^ = — y /f | ^~ L ^ 1 ~ L ) . Later in this subsubsection, we analyze the explicit behaviors of the 
scalar field and the Hubble parameter for the cases with 4> ~ and <f> — > —00 successively. 

We first consider the behaviors of the solutions around cf> ~ <f> c . Note that G(4>) is positive definite. Therefore, the 
velocity <fi goes to zero only when L(cj)) = where <fi — <p c . In the past, the scalar field monotonically increases with 
time until it reaches <p c . Thereafter, it monotonically decreases. As t — >• ±00, the scalar field goes to 4> — >• —00. The 
Hubble parameter at <fi = <j) c becomes 

Near L ~ 0, the Left-Hand-Side(LHS) of Eq. (JT3J| is expanded to be 1 + ^-L 2 + ■■■ and gives L ~ 6 1 / 4 ^/3(^) c - <t>). 
Integrating Eq. (|15[) around <fi c , we get 

where we assume that 4> = <fi c at t = 0. Because the scalar field behaves monotonically before and after (f> c , we may 
use the field value 4> as a time as follows: 

i-4> c , t<0, 
> c -</), t>0. 

In terms of this time, the metric can be written as, 



ds 2 = - 2G 2(^ L 2(0) ^ 2 + a\<t>){dx 2 + dy 2 + dz 2 ), 



2 The parameter f3 in this work corresponds to A/ y/6 in Ref. 1 1 -tH . 
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where the scale factor is 

a(4>) = doe 75 



bit' 



-<!L 



(l-L) 1 +' 9 (l + i) 1 -f 



2x/3(l-/3 2 ) 



We next consider the asymptotic limit t — > ±00 with <fr — > —00. The asymptotic behaviors of the scalar field and 
the scale factor are dependent on the value of (3. In the case of < (3 < 1 one can show that the solutions asymptotic 
to the fixed point solution (fTTj). The procedure is as follows: Noting the fact that the Right-Hand-Side (RHS) of 
Eq. (TT3")) is very small, one sets L ~ /3(1 — e). Then, from Eq. ([TB")) . one gets 



(l-/3)^-(l + ^)V eX p 



31-/? 2 

/3 



(0 - 4>c) 



Then, one can integrate Eqs. (TT2")) and (fTSj) to get the evolutions of the scale factor and the scalar field asymptotic to 
Eq. (HB). 

For (3 > 1, the RHS of Eq. (TT3")) becomes very large for — > —00. To achieve this behavior, we may have L — >• 1 
or L — > — 1. First consider the case with L ~ 1. In this case (/> < and corresponds to the later time universe with 
t — > 00. Explicitly we have 



L -} 1 



2ff 



3-1 



As £ — > 00, the scalar field and the scale factor behaves as 



<t>(t) -> (/3-l)0 c - 

tt(t) -> ^(t + ty) 1 / 3 



^log 



22(3 + 1) 



(1-/3-1)^ 



y|log(VA(t + */)), 



(16) 



The asymptotic value explicitly dependent on the choice of <p c and deviates from the fixed point solution (ITT1) . Wc 
next consider the case with L — > — 1. Explicitly we have 



L -> - 



1 2 ^ 

(1 + /3- vs(g + i)^-^ e ) 
~ 3+1 c 



23^ 



In this case, </> > and corresponds to the early universe. As t — > —00, the scalar field and the scale factor behaves as 



4>{t) ->■ C3 + i)0 c + W-iog 

o(t) -> a (t-t ) 1/3 , 

where £ denotes the initial time of the universe. 
The Universe expands with accelerating rate if 



3+1 

2 2 (3-i) 



y|]og(VA(t-to)), 



(17) 



1-3L 2 
3(1 -i 2 ) 



> 0. 



which leads to |L| < l/y/3. Therefore, if the universe starts with a small L (< 1/3), then it will experience temporal 
inflation. 



IV. EXACT SOLUTIONS AND THEIR PROPERTIES 



In the cases of the potentials other than the constant and exponential, wc fail to find the whole solution spectra but 
obtain a few specific solutions for each potential. For a given equation of motions, it is important to understand the 
role of fixed points of generating equation (|6]) on the evolution of the universe and the stability of the specific solutions 
which may or may not asymptotic to the fixed point. In this section, we obtain the condition that a generating function 
is asymptotic to fixed points. We next find the specific solutions corresponding to the generating function for the 
cases of the power-like and double-well potentials. 
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A. Stability of generating functions 



For a given scalar potential, the equation of motion ([3]) and (j4]) has a fixed point </> c , which can be obtained from 
equation G'(4> c ) = i.e., G(<fi c ) = y/V{<j) c ). To see the properties of the fixed point, we study perturbations of the 
generating equation (j6)). Let Go(<p) satisfies Eq. (j6)) and is asymptotic 3 to the fixed point for the potential V (</)). To 
study the stability we perturb the generating function 

G(0 = G o (0(l + e(0). 
From Eq. ([5]), the perturbation e should satisfy, to the first order, 



<4>) 

Therefore, the perturbation takes the form, 



e(4>) = e exp 



2G O (0) GoW 



S((j)')d<j)' 



(18) 



Now, we find that the perturbation goes to zero at 4> c iff J S(tf>')d4>' 
the function S(<j>) behaves as 



S(<t>) 



-4>c) 



-co. This behavior can be accomplished if 



with n > 1, s 2 > 0. 



(19) 



On the other hand, if s 2 < and n > 1, e(</>) becomes divergent at <f> c . In this case, <f> c corresponds to an unstable 
node, which we are not interested in this work. Assuming s 2 > and n > 1, there are two ways for the generating 
function to satisfy Eq. (fTJ 



1) 



Go 2s 2 (0-0 c ) 



2): ^° 
J Go 



G 3 10-^1"+!' Go |(/»-0 c |"+ 1 

Integrating the equation for each case, we get the behavior of the generating function Go (<f>), 

"3|0-0 C |" +1 " 



(20) 



1) : GaW 

2) : G o (0) 



s/A exp 
•\/Aexp 



2 (n + l)s 2 

o2 



(n-l)|0-^ o | T 



G o (0) = VA|^- c p s with 



(21) 



The case 1) corresponds to a solution approaching to the fixed point. We deal this case in detail in Subsec. IIVBI For 
the case 2), the generating function Go{4>) goes to infinity as <p — > <f> c , which implies the Hubble parameter diverges at 
the field value. In addition, the kinetic energy of the scalar field is singular since G' Q (<p) — > oo as <fi —> <j> c . Therefore, 
the space-time will develop a singularity and should end (or begin) there. In addition, a well-defined scalar potential 
may not diverge for a finite scalar field value. Therefore, we will not discuss the case 2) any more. 



B. Power-like potentials 

As announced in the previous subsection, we consider the first case of Eq. ([2~T]) . Without loss of generality, we may 
set 4> c — 0. Since we are observing the behavior of the generating functions around <f> c = 0, we may use 

G(0) = v / A(l + ^-|C +1 ) (22) 



3 The scalar field and the evolution of the universe has a fixed point there. However, in this work, we mention that the generating function 
is asymptotic to the fixed point for simplicity. 
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prefer to the first equation of Eq. (|2T|) . For <fi = to be a stable fixed point, the potential needs to satisfy « = j^y > 
0, n > 1. In this subsection, we do not restrict the values of u and n but discuss the general properties of the 
solution generated from (|22j) . Even though we use the absolute value in the equation, we restrict the motion of the 
scalar field to the region <j> > because the velocity 4> goes to zero there for most cases we are interested in. 
The scalar potential obtained from the generating function (|22j) is given by 



V{4>) = A(l + -^-|0r +1 ) 2 - |a m V". (23) 
n + 1 3 

For n < 0, the potential goes to negative infinity as </>—»• 0. For n = 0, after the change of variable (j>+ (n + 1)/" — > 0, 
the potential becomes nothing but a 'mass term+constant' with its mass-squared, m 2 = 2A/x 2 , and the constant part, 
Vo = —^m 2 . For < n < 1, the potential has a local maxima V(0) = A>Oat0 = O and two degenerated minimum 
at ±4> m satisfying (1 + u</>^ +1 /(n + l))(f>]^ n = For n = 1, the potential describes nothing but a A</> 4 theory, which 
will be dealt in Subscc. riVCl in detail. For n > 1, it has several local minima and maxima at (j> = and at the values 
of ±0 m . The first derivative V'{4> = 0) is zero, singular, or finite for n > 1/2, n < 1/2, or n = 1/2, respectively. 

Let us discuss the general properties of the solution starting from the case with fi < 0. Recalling H cx G(</>) and 
cx —G' (</>), the universe will be in a contracting phase for > (— /_t/(n + l)) 1 /™. Because <fi > at all time for 
7i > — 1, the expanding universe at <fi = enters into contracting phase for > (— [i/(n + l)) 1 /". This behavior 

of the universe can be understood by noting that the equation of state w = — 1 + f (i+ M /(n+i)x0' i + 1 )^ diverges as 
0^ (-M/(?i + l)) 1/(n+1) . 

In the rest of this work, we assume that the Universe expands forever (/i > 0). The acceleration of the scale factor 
becomes 



M ( I _ V2| 0r + ( 1 + V5|0|- + ^tI^I" +1 ) 

3\/j n+1 n + 1 / 



The Universe expands with accelerating rates always except for a possible short decelerating period which exists if 
a > n+1 

Without loss of generality, we confine <fi > since the dynamics will end (or at least bounce back for n > 0) at 
(f> = 0. Then, the scalar field and the scale factor for n ^ 1 evolve as 



a(t) = a cxp I Hit - 2(n+1) (2(n - l)A*ff/t) "- 1 J • (24) 

The solution for the n = 1 case is given in the next subsection. 

For n < 1, the domain of time is (— oo, 0). It takes a finite time to arrive at <f> = 0. For n < 1/2, the space-time will 
develop a singularity at t = because V'(0) diverges. For 1/2 < n < 1, noting the behavior of V'(0), the space-time 
will be regular and can be extended to the region with positive t. There, the scalar field bounces to increase and scale 
factor expands, however, we can not exactly describe the later evolution with the present formalism. For n > 1, the 
domain of time is (0, oo). The scalar field asymptotically approaches to 4> = and the Hubble parameter approaches 
to H] leading the Universe to eternally inflating phase. 



C. A<?!) 4 theory 

Now, let us examine the A</> theory with n = 1 case in the previous subsection. The scalar potential is given by 



V&) = \(tf- <P m 



Vn 



where, with Hi — and A = A/i 2 , 



/ 2 Hi \ _ 2A / Hi 



(25) 



Note that the three parameters A, (f> m , and Vo are not completely independent from each other but are constrained 
by the last equation of Eq. (f2~5"T) . For Hi < |^/A/3 (c/) 2 n > 0), the potential becomes a kind of double- well potential 
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with its minima Vq located at <p — ±fa n . For Hj > |a/A/3 (0 2 t < 0), it becomes a single- well potential with a unique 
minimum A at 4> = 0. 

The scalar field and the scale factor are given by 

fat) = fae~ 2 ^\ a (t) = a cxp Ujt ^e'^^ . (26) 

The domain of time runs (— oo, oo). In the past, the value of the scalar field is decreasing from infinity. The scale factor 
of the Universe is expanding with accelerating rates which may not be related with the inflation we are interested in. 
As time increases, the scalar field continually decreases and approaches to zero in the future. The scale factor will 

experience a decelerating period during \/2(l — ij 1 — Hi/y/X/3) < 4> < \/2(l + \J 1 — Hi/^X/3) if Hj < y/X/3. 
Eventually, the Universe goes into eternally inflating phase as the scalar field approaches to the fixed point = 
with its Hubble parameter Hj. 




FIG. 1: The evolution of the scalar field with perturbation. The scalar field starts to evolve from the initial value 
arrive at <j>f. Later, it will rolls down over the potential and will oscillate around cj) m . 



and 



As mentioned in subsec flV A[ the generating function G(<p) is asymptotic to a fixed point of the Einstein equation. 

Therefore, small deviation from G{fa) = VSHi + ^-fa 2 tends to decrease as 4> 0. Once the scalar field arrives at the 
fixed point, the Universe is destined to inflate eternally driven by the scalar potential at cf> = 0. Note however that 
= can be an unstable equilibrium if 4% n < 0. In this case, any small perturbations will make the scalar field roll 
down and the inflation stop, which expectation appears to contradict the stability analysis in Eq. (fTT)|) . In fact, the 
stability (|19[) does not directly imply that the eternal inflation is a later time attractor in the solution space because 
the generating equation (jB]) is a differential equation with respect not to the time but to the scalar field. Now let us 

check how the stability affects on the time evolution of the scalar field. From Eq. (f2"2"j) , with Gq{4>) = + -^0 2 , 

we get the small perturbation (fT8|) around Gq(4>) evolves as 



<4>) 



2\/a73 2Hj 



2Ht 



■ cxp 



Here fa (> fan) is the field value at the moment </> = as in Fig. Q] and is determined by the condition, 



= G'{fa) = G' (fa) 



3e G o (0) 2 
2 G' (fa 



Vxcj) 



3e (2Hi + yfX/Zfa 2 
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(A/3)02 



(27) 



At 



it determines the value of eg to be 



(A/3)0 2 



3 (2Hj + v/V30 2 



o(fa 
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One may seek for another solution of = around 0^0. If e(0) decreases faster than 1 as — ^ 0, there does not 
exist root of Eq. (|28|l other than = and this is the case if 

2 < 0. 

Note, in this case, that the potential has a minimum value at = and the fixed point solution now plays the role 
of a later time attractor. 

In the case of a double- well potential with 4> 2 m > 0, there exist another solution 0/ of Eq. (|2T|) . satisfying 




e -sf/ 



2-/V3 -2,4? 



2Hi 



(28) 



Assuming 0j ^> 1, 0f is given by 



4> f ~ 0i exp 



30? 



16(1 



3Hj 




(29) 



where the value 0/ is suppressed by both of the exponential factor and the inverse power. In conclusion, the per- 
turbation on the initial condition, where the scalar field starts to evolve at </>,;, makes the scalar field fail to evolve 
to — but to arrive at 0/ which value is extremely small. Therefore, the Universe fails to go into the eternal 
inflating phase but still experiences a very long period of inflation when the scalar field stays around / . Note that 



the accelerating expansion starts when the field takes the value 0i = V2~(l — y 1 — Hj / \/A/3 ) . After the time when 
the scalar field reaches to 0/, the field value will rolls down over the potential. During this period, the universe still 
inflating until the scalar field value reaches the prescribed value. After that, the universe enters into the reheating 
regime. Therefore, the total number of e-folding during the accelerating expansion of the universe before it enters 
into the reheating phase is roughly, 



iV 



2H I (t f -t 1 ) 



log 



<Pl 



30? 



3/4 



(30) 



where we have ignored subdominant contributions. The long period of inflation is also noted to happen for the 
inflection point scenario where the potential locally has a cubic scaling [26| . 



D. Double well potential with hyperbolic function 

It is interesting that the appearance of a long inflating period around the fixed point of Eq. © even the scalar 
potential is in a broken phase around the fixed point. Since the generating function (|22p asymptotic to a fixed point 
comes from a general argument in Eq. (|19|1 . the long inflation appears to happen quite generally in the evolution of 
the universe including a scalar field. In this subsection, we present another example. Let us consider the hyperbolic 
generating function, 



G(0) = VA(cosha0 - 0). (31) 
Since we are considering an expanding Universe around = 0, we require /3 < 1. Now, the scalar potential becomes 



V{4>) = A 1 



2o2 
3 



cosh aQ 



V , 



2a 2 A 



\ x l-2a 



satisfying cosha0 m = — 



For the potential be bounded below we assume a 2 < 3/2. The potential has a 

■ - l) 2 at = 
and if B > 1 



where V - — ^ - 1 _ 2a!i/3 / 

unique minimum with V(0) = A(/3 — l) 2 at0 = Oif/3<l — |a 2 and has two degenerated minima V(0 m ) = Vq at 

P 



la 2 
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The scalar field and the scale factor behave as 

(j)(t) = — logcoth (a 2 -ffji) , 

a 

a(t) = a e- fmit [sinh (2a 2 i?/i)] ^ . (32) 

where we assume that the scalar field rolls down the potential from the positive side initially. The domain of time is 
(0,oo). As t — > 00, the scale factor exponentially increases with its Hubble parameter H(t) — > (1 — f3)Hj. At both 
ends of the time, the Universe expands with accelerating rates if B < 1. An intermediate decelerating period exists if 
B 2 + 2a 2 > 1. 

Similarly as the case of Xcj> 4 theory, the generating function (f3Tj) approaches to the stable node at <f> = 0. Let us 
check whether the corresponding solution (|32|) play the role of later time attractor or not. From Eq. ([3~Tj) . we get the 
small perturbation (|18j) around the generating function is given by 



e(</>) = £0 



cosh afa — 8 
cosha0 — (3 



sinh a(j) ( tanh 



sinh afa \ tanh 




Here we choose </>i (> </> m ) to be the field value at the moment = 0, which is determined by the condition 



= G'U) = G'U) 



3e G (4>f 



2 G' {<t>) 



%/Aa sinh ac 



3e (coshatfi — 0) 
2 a 2 sinh 2 ad> 



(33) 



At 



it determines the value of eo to be 



2a 2 



sinh 2 1 



3 (cosha^i — /3) 2 



Note that the value eg will be of order a 2 . Because we are doing a perturbative calculation, we simply assume a 2 « 1. 



Now, let us seek another solution of <fi = around 4> ~ 0. If e(</>) decreases not slower than 
not exist non-zero root of <p = other than <p = satisfying cf> < <p m and this is the case if 



0. there docs 



l-/3> 



2a 2 



which corresponds to a single- well potential. For the case 1 — j3 < , Eq. p3p is given by 

(sinh acf> f ) 2 ~ 5^ (tanh ^ ) ^ (sinh f" & (tanh ^ ) ^ 



(cosha^/ — /3) 



Assuming a<f>f -C 1 and a</>,; S> 1, we get 



(cosha^i — (3) 



acj)f 



8 3(1-/3) 



exp 



3 -2a 2 



4a 2 -3(1-/3) 



which is exponentially suppressed. Since we are assuming a -c 1, this equation holds for /3 ~ 1. 

Therefore, the scalar field will evolves to 4> = with time and the Universe goes into eternally inflating phase in 
later times even in the presence of initial perturbations. 

In summary, the eternal inflation is a later time attractor if the potential has a local minimum at <p = for various 
potentials including the X(j) 4 and hyperbolic potentials given from the generating functions having the behavior (|22[) . 
On the other hand, if <f> = is in a broken phase, it cannot be a later time attractor not because the perturbation 
does not vanish there but because the field fails to arrive at <j) = 0. However, the property of the generating function 
as asymptotic to the fixed point leaves imprint: the Universe experiences a very long inflating period, during when 
the scalar field stays around its minimum value </>/ ps 0. In this model of inflation, we do not need to impose artificial 
conditions such as the slow-rolling which constrain the motion of scalar field to give enough e-folding for the inflation. 
The long period of inflation is nothing but a consequence of the evolution of the preinflationary Universe. 
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E. 



potential 



Let us consider the potential with </> 6 order interactions. We display this solution as an example of exact solution 



where a scalar field settles down at a stable minima of the potential at 
considering, is 



G{4>) 



1 + 



0. The generating function, we are 

(34) 



The potential is given by 



V{ct>) 



A ( 1 + ^ 2 



— 2 (/i + a0) 2 . 



The change a — > —a is the same as <p — > —<j). Therefore, we may restrict a be positive definite without loss of 
generality. The derivative of the potential becomes 



V'(4>) = 2A(t)(n + aq 



3 3^ 2^ 



Therefore, <j> = and <j> = —fi/a are the place where the derivatives of the potential vanishes. Around <f> = 0, 
V(</>) ~ A(l + /i(l — 2/i/3)0 2 + •••)• Therefore, /i = is a stable equilibrium or an unstable local maximum, 
respectively for < /i < 3/2 or /i > 3/2. For \i = 3/2, <j> = becomes a saddle point of the potential satisfying 

4 + a) 2 H • Therefore, (j) = -fi/a 



A(l + ^)-M 1 + f + 



V'(0) = = V"(0). Around = -/z/a, V{4>) 

is an unstable local maximum. 

Noting <f> = —2Hj(f>(fj. + a<f>), we find that the point <p = —fJ,/a is an unstable node and 
Now, the evolution of the scalar field and the Hubble parameter becomes 



is a stable node. 



4>{t) 

a(t) 



Ae 2H '^ -a 



a,Q exp 



(Hi 



6a 2 



\2a 2 Hj \Ae 2H '^ -a (Ae 2H '^ - a) 2 



log \Ae 



2H!fj,t 



(35) 



where A is an integration constant. If the scalar field starts to evolve from the value 4n > —fi/a, it will arrive at = 
asymptotically, where the unverse will expand exponentially. On the other hand, if it starts from values </>,; < —/j,/a, 
it evolves to negative infinity. A 6 potential was also studied in Ref. [2(| with the fake supergravity method, in which 
they show that the fixed point can be a local maximum of a minimum of the potential considering a phantom scalar 
field. The ways were extended to include models with two scalar fields in Ref. (2lj . 



STABILITY ON THE PERTURBATIONS OF SCALAR POTENTIAL 



In the previous section, we show that many scalar potentials derived from the generating function of the form 
leads the Universe to experience a long inflating period. It is interesting how the solution changes if the potential 
does not take the form derived from Eq. (|22|) . To examine this possibility, we perturb the potential to get a new 
solution around the known one. Let Go(<p) be a generating function of the form (f2"2")l for a given potential V(4>). We 
try to find a new generating function which may or may not takes the form (|22p and corresponds to a new potential 
V(<j)) + v(<j)) with v(4>) *C V((j>). New solutions of the generating equation nearby Go(<^>) can be found by setting 

G(4>)=G (4>)+e(<t>). 

Equating to first order from Eq. ©, e satisfies Go(4>)s((f>) — ^G' (<f))e' (4>) = -^p, which can be integrated to give 



e(<p) = Q{ct>) 



Cl 



Q(4>) = exp 



aW) d<t> ' 



where c\ is an integration constant determined from an initial condition, which is related to the change of the initial 
condition at initial time. If we are interested in the modification of the solution due to the potential change, we may 
set Ci = as in the rest of this section. 
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Now, as an example, we examine the Xcj) 4 potential case by adding a small constant potential term v(<j>) = v. Then, 
we have 



e{(j>) = e 8 



3^2 3u 3 2 

— i 3H i 
8VA 1+ 7JT, 



where E a (y) denotes the exponential integral function. Noting the series expansion of the exponential integral function 
around y = 0, E a (y) = y a T(-a) + ± + ^ + + 0(y 3 ), we find that 



e(t/>) = 



3t; 
8^ 



_3Hi 



3Hj_\3Hi 



o(</> 4 



(36) 



Including e(</>), the generating function becomes 
G{<j>) = V3Hj 1 1 



v \ A 

2^/(4^/173 + 3i?/) i 2 ( 



which has of the same form as ([22)1 with n = 1. Therefore, with this perturbation of the constant addition of the scalar 
potential, the = still plays the role of a fixed point and the Universe will inflate for a long time with modified 
Hubble parameter when the scalar field arrive around = 0. 



VI. SUMMARY AND DISCUSSIONS 



We presented an algebraic method to find solutions of the Einstein-scalar field equations in spatially flat Fricdmann- 
Robcrtson- Walker space-time. The equation for the scalar field part turns out to be integrable by choosing the Hubble 
parameter as an appropriate function of the scalar field and its time-derivative. As a result, we find that the Einstein- 
scalar field equation is equivalent to a generating equation: a nonlinear differential equation (HU) which relates a 
generating function G(cj>) with the scalar potential. For the equation, we have displayed two different interpretations 
as a differential equation to find the generating function and as an algebraic equation presenting the scalar potential 
from the generating function. 

Following the first interpretation, the solution spectra were obtained in the cases of constant and exponential 
potentials, which reproduces the results in Rcf. [22J. However, the differential equation ([5]) is hard to solve in 
general, which makes the second interpretation powerful. Following the second, we find several exact solutions of 
the Einstein-scalar field equation for various scalar potentials such as the Xcj> , the power law, and the double- well 
hyperbolic functions. We can control the behavior of the solutions since the derivative of the scalar field and the 
Hubble parameter are simply proportional to G'(cf>) and G(</>), however, we loss the power to control the potential 
as a return. We have analyzed the stability of the generating equation ^ and found that a fixed point exists if the 
generating function takes the form around </> c 

G(0)~ v/A(l + ^-|0- C |" +1 + •••); »>1, M>0. 
n + 1 

The simplest n = 1 case provides the A</> 4 theory with its zero energy is determined by the coupling and the mass 
term. 

The solution corresponding to above generating function were presented. If the mass-squared is non-negative at 
the fixed point, it was shown that the fixed point becomes a later time attractor of the solution. On the other hand, 
if the mass-squared is negative definite (double- well) , the scalar field fails to arrive at the fixed point but bounces 
back just before the fixed point. Even though the fixed point is not a later time attractor, there exists a very long 
period of inflation during the period when the scalar field stays around the fixed point. We additionally examined 
the possibility whether the small change of the potential may alter the later time property or may not. We found 
that a constant addition on the scalar potential won't change the property of the generating function. Therefore, the 
Universe will experience a long inflationary period even if the potential does not take the form generated from the 
above form but deviates a bit. This proves that the appearance of a long inflation through the evolution of a scalar 
field on top of a local maxima of scalar potential is a frequent occurrence in the Universe. 
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